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ABSTRACT 

We use the well tested StarTrack binary population synthesis code to examine the properties 
of the population of compact object binaries. We calculate the distribution of masses and 
mass ratios taking into account weights introduced by observability in gravitational waves 
during inspiral. We find that in the observability weighted distribution of double neutron star 
binaries there are two peaks: one for nearly equal mass systems, and one for the systems 
consisting of a low and a high mass neutron star, q = 0.6 — 0.7. The observability weighted 
distribution of black hole neutron star binaries is concentrated on the systems with the mass 
ratio q — 0.3 — 0.5, while for the double black hole binaries the observability weighted 
distribution is dominated by the massive, nearly equal mass binaries with q > 0.7. 
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1 INTRODUCTION 

We are currently witnessing a large incr ease in the sensitivity of t he 
gravitational wave observatories. LIGO iAbramovid_£talj l99jHs 
alread y taking data, the development of VIRGO I Bradaschia^ydJ 
|199(J) shows gre at advances, GE O600 iPanzmann & et alJI 1992ft 
and TAMA300 (Tsubono 1995) are operational. In the coming 
years an even more sensitive Advanced LIGO will begin taking 
data. Out of a number of potential sources of gravitational radia- 
tion the most promising are probably mergers of compact object 
binaries, i.e. binaries consisting of black holes (BH) and/or neutron 
stars (NS). These are the only sources for which observations in 
the electromagnetic domain are consistent with emission of grav- 
itational waves. Present efforts to examine data from gravitational 
wave detectors show that such detections rely heavily on availabil- 
ity of accurate templates. This provides a case for the importance 
of accurate merger calculations. The data analysis relies on cross 
correlating the data with a number of templates. Scanning large 
volume of the parameter space requires using a large number of 
templates and may hinder detection of a real but low amplitude sig- 
nal. Any possibility to limit the amount of templates required or to 
show in which region of the parameter space a detection is most 
likely may improve chances of seeing the gravitational waves. 

Thus it is important to ask the following questions: what are 
the most likely objects to be observed? what are the most impor- 
tant parameter sets to explore? In order to answer them one needs 
to investigate the properties of population of compact object bi- 
naries. Observations provide us with six do uble neutron star bi- 
naries fThorsett & Chakrabar tvlll999l:lBurgav et alJl2003l) . The ra- 



dio selected sample of double neutron star binaries is biased to- 
ward long lived systems. However, we do not know any black hole 
neutron star nor double black hole binary. Therefore inferring the 
properties of the population of compact object binaries solely on 
the observations of these few systems may lead to erroneous re- 
sults. A different approach - the binary population synthesis al- 
lows to investigate the properties of such systems from a theo- 
retical point of view. Binary population synthesis requires how- 
ever, a thorough investigation of the systematic uncertainties due 
to parametrization of various stages of stellar evolution. The pop- 
ulation synthesis studies have already been used to estimate the 
rates and properties of the merg ers that can be observed by the 
gravitational wave obse rvatories iLipunov et al. 1997b; Fryer^UilJ 
19981 IPorteeies Zwart & Yungelsonlll998l: iBethe & Brownl | l998|: 
Fryer et alj Il999t iBelczvnski & Bulikl Il999t iBulik et alJ fl999 



Belczvnski et alJl2002d: iNutzman et al1l20oj) . It has been shown 



that the observed sample will most likely be dominated by the 
mergers of double black hole binaries iLipunov et alJ Il997at 
Buli k & Belczvnskil 120031) The distribution of observed chirp 
masses was found to be a very sensitive indicator of the stellar 
evolution model while being relatively not se nsitive to the star for- 
mation rate history and cosmological model iBulik et all2004l) . A 
preliminary study of the distribution of m ass rati os in compact ob- 
ject binaries was presented by B ulik et all 1200 3l) . 

In this paper we use the StarTrack population synthesis code 
to investigate the distribution of masses and mass ratios in the pop- 
ulation of compact object binaries. We use a convention where the 
mass ratio q in a binary system is defined as the ratio of the lower 
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Table 1. Population synthesis models. We list the number of coalescing compact object binaries produced in each simulation. For detailed models description 
see § 2.1 and § 2.2. 



Model 


Description 




N produced 


A 


standard model described in § 2. 1 




5761 


Bl 


zero kicks 




21535 


B7 


single Maxwellian with c — 50kms ^ 




17747 


Bl 1 


single Maxwellian with cr ^ 500, kms ^ 




2155 


B13 


Paczvnski 1 1990) kicks with Vk = 600km s~ 1 




8270 


c 


no hyper— critical accretion onto NS/BH in CEs 




4798 


El 


CE efficiency' ckqjh] x A — 0.1 




894 


E2 


CE efficiency! ct^cE x A = 0.5 




3489 


E3 


CE efficiency: ckce X A = 2 




8504 


Fl 


mass fraction accreted in non-cons. MT: f a = 0.1 




2483 


F2 


mass fraction accreted in non-cons. MT: f a = 1 




4644 


Gl 


wind decreased by / w i n( j =0.5 




9395 


G2 


wind changed by / wind = 2 




5517 


J 


primary mass: <x M-j -2 ' 35 




8220 


LI 


angular momentum of material lost in non-cons. MT: j 


= 0.5 


6660 


L2 


angular momentum of material lost in non-cons. MT: j 


= 2.0 


5547 


Ml 


initial mass ratio distribution: &(q) <x q~ 2,7 




852 


M2 


initial mass ratio distribution: <£>(<?) <x <j 3 




11225 


O 


partial fall back for 5.0 < M C < 14.0 Mg 




4116 


s 


all systems formed in circular orbits 




4667 


Zl 


metallicity: Z = 0.01 




5199 


Z2 


metallicity: Z = 0.0001 




7074 



mass component to the higher mass one and therefore is always less 
than unity. In section 2 we shortly describe the code, and demon- 
strate the difference between the volume limited and flux limited 
distributions of masses of compact object binaries. We present the 
results in section 3, and conclusions in section 4. 



2 CALCULATIONS 

We are using the StarTrack b in ary ev olution code described 
in detail by Belc zynski et alj l2002d) . The code is well 
tested and has been used in various astrophysical applica- 
tions: analysis of gamma ray-burst progenitors I Belczvnski et al. 
2002b), tracing of evolutionary history of individual bina- 
ries iBelczvnski & Bulildf2002ll. investigation of the mass spec- 
tra of compact objects l Belczvns ki et alj 1200231. StarTrack 
population synthesis code was specifically designed to cal- 
culate the merger rates and physical properties of com- 
pact object binaries. It was compared with a several other 
codes fe^jLijDuno^e^aij j^ 

<1998l) Ide Ponder & Vanbeverenl il998h . iFrveretalJ J 19991) . 
Nele manset alj fcOOlhT . The comparisons showed some differ- 
ences, however they were understood within the different model 
assumptions. Since StarTrack was designed to deal mostly with 
systems containing NSs and BHs, our input physics was updated 
and revised as compared to the other codes in respect to compact 
object formation. As a result we have recognized new evolution- 
ary NS-NS formation scenarios, and we have shown that massive 
stellar BH may dominate the population of double compact objects 
(~ 10A/q) observed in gravitational waves. Last, but not least, our 
predic tions of NS-NS Galactic coalescence rates iBelczvnski et alj 
2002c) are in good agreement with the most recent constrai nts ob- 
tained from the observed sample of these systems I Kalog era et alj 
l200ltlKim et all2003l) . 



2.1 Standard Model 

Within the code the evolution of single stars is parametrized by the 
modified formulae of Hur ley et alj 1200(3) . The single star evolu- 
tion includes such stages as the main sequence, evolution on the 
Hertzsprung gap, red giant branch, core helium burning, asymp- 
totic giant branch, and evolution of helium stars. Two major mod- 
ifications include low mass helium star evolution and calculation 
of compact obj ect masses. In particular, following a number of 
studies ff)elgado & Thomasll98ltlHabetsil987lJAvila-Reesell993t 
Woosl ev et alll9 95) we allow the low mass helium stars 4Mq) 
to develop deep convective envelopes. Presence of convective en- 
velope plays an important role in the behavior of the donor star 
in the Roche lobe overflow event, and may eventually lead to the 
development of dynamical instability and common envelope (CE) 
evolution, and possible tightening of the binary orbit. 

The original Hu rley et alj yOOO) formulae are used to calcu- 
late the final CO core mass of a given compact object proge nitor at 
the time of supernova/core collapse event. We use IWooslevU l986) 
stellar models to obtain mass of the final FeNi core correspond- 
ing to a given CO core mass. The FeNi core is collapsed to form 
proto-neutron star, and then we use the results of core collapse hy- 
drodynamical calculation of Fryer (1999) to calculate the amount 
of fall back material and the final mass of the newly formed com- 
pact object. We use the following algorithm to derive the masses of 
a newly formed compact object M rem : 

C Mfoni Mco < 5 M 

M rcm = t M F eNi + fib(M - M F eNi) 5 < M co < 7.6 (1) 
{ M Mco > 7.6 M 

where Mf ni is the mass of the FeNi core, Mco is the mass of the 
CO core, M is the total mass of the star prior to the explosion, and 
/fb is the fall back factor, < /ft, < 1 depending on the mass 
of the star. This simple formula represents well the results of de- 
tailed numerical calculations. We do verify the sensitivity of our 
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results to changes in the particular numerical values in equation^ 
see e.g. model O below. Varying stellar evolution parameters like 
the strength of winds, or metallicity leads to different core masses 
for a star of given initial mass and also alters the initial final mass 
relation for single stars. We find that NSs are formed without a sig- 
nificant amount of fall back material, while BHs are formed either 
directly (prompt collapse of a massive star) or through partial fall- 
back of material onto proto-neutron star. 

The binary evolution takes into account orbit changes due to 
wind mass lo ss, and tidal interactions. Wind mass loss rates are 
adopted from Hurlev et al. 1 2000 ) and they depend on the stellar 
parameters of mass losing component (its composition, mass and 
evolutionary stage). Specific mass loss rates are adopted for naked 
helium stars, the luminous blue variables and pulsating stars. In sta- 
ble mass transfer (MT) calculations we allow for non-conservative 
evolution. We assume that part (f a ) of the transferred material to 
accreted onto companion star, while rest is ejected from the system 
with the specific angular momentum (J, expressed in the units of 
the binary angular momentum). In our standard model we adopt: 
f a ~ 0.5 and j — 1. If the RLOF episode is dynamically unsta- 
ble, we follow the spiral in through the common envelope phase. If 
system avoids the merger, we calculate the final orbital separation 
using standard energy conservation based prescription of IWebbinkl 
1984). The evolution through CE phase depends crucially on the 
efficiency of the orbital energy input into the donor envelope (a ce ) 
and the specific binding energy of the envelope (A). Only the prod- 
uct of these two largely uncertain quantities enter the calculation, 
and we use q co X A = 1 in the standard model, however we also 
check the sensitivity of the results to this parameter. During the CE 
spiral i n we al l ow for hyper crit i cal ac c retion onto NSs a nd BHs 
fe.g..jBlondirJ <1986MChevalie3 ll989h . IChevalieJ <1993l)lBrownl 
1 1995)). As a result several tenths of solar mass may be accreted 
onto the compact object, and in particular the top heavy NSs may 
collapse and form BHs. The full description of the hyper c ritical 
accretio n treatment is given in the appendix of Belczvnski et al. 
<2002d) . 

Supernovae explosions are treated in detail. The explosion 
takes place at the randomly selected place on the orbit. We allow 
for explosions on the eccentric orbits, for uncircularized systems. 
We take into account the instantaneous mass and angular momen- 
tum loss form the binary system. Also a natal kick is added to 
the orbital velocity of the newly born compact object to account 
for the SN asymmetr y. Kicks ar e sele cted from the bimodal dis- 
tributions of ICordes & Chernofj <1998h . a weighted sum of two 
Maxwellians, one with a = 175 km s _1 (80%) and the second 
with a — 700 km s _1 (20%). A binary is either disrupted in the 
explosions, in which case we stop the evolution, or if it survives we 
follow the evolution on the new binary orbit. 

The initial mass of the primary M\ ams is drawn f rom a power 
law initial mass function (IMF) distribution oc M~ ~ 2 - 7 <Scaloll986h 
within the range 8-100 Mq. The secondary mass is obtained as 
Mz ams = qM\ a rnx , where q is the mass ratio and is drawn form 
a flat distribution jKuipeJl935l) . We allow for eccentric initial or- 
bits, and the eccentricities are drawn from a thermal distribution 
oc e lHeggiell975l : lDuauennov & MavoJl99ll) . Fi nally, the o rbital 
separation distribution is taken to be flat in log a (Abt 1983), and 
separation are chosen from few (so the stars are not formed at the 
contact configuration) up to maximum of 10 5 solar radii. We evolve 
our stars for a maximum Tnubbie = 15 Gyrs. The evolutionary 
model described above is chosen as our reference (standard) model 
and marked with latter 'A" on the following figures and tables. 



2.2 Parameter Study 

In order to asses robustness of the results we investigate 20 extra 
different models of stellar evolution, where we vary the parame- 
ters describing various stages of stellar and binary evolution. The 
models used are listed in Table 1 . The range of models represents 
the current state of knowledge and uncertainties about the binary 
evolution. All models are calculated with 2 x 10 6 initial binaries 
each. 

In models marked with letter "B" we vary the distribution 
of natal kicks compact objects receive when they are formed . 
This is rather uncertain part of evolutionary model as we still do 
not kn ow the mec hanism behind the SN/core collapse asymmetry 
(Buras et al. 2003) . Therefore, we change the kicks quite drasti- 
cally, from rather non-realistic model with no kicks (B 1) to the very 
high kicks of model Bll. The higher the kicks, the less compact 
object binaries we form, since the higher kicks tend to disrupt the 
progenitor systems. This is one of the most important parameter as 
far as the number of compact object binaries are concerned (close 
to an order of magnitude change). 

Since the CE evolution is another highly uncertain part of our 
evolutionary scheme, in models "E" we change the efficiency with 
which orbital energy is transformed into the unbinding the enve- 
lope of the donor star, while in model "C" we turn off the accretion 
onto compact objects during that very short lived phase. In mod- 
els with small CE efficiency (El -2) it is found that the number of 
compact object binaries are significantly reduced. This is due to the 
fact, that many binaries, evolving through the CE phase, will merge, 
thus aborting compact object binary formation. On the other hand, 
increase of the efficiency (E3) or shutting down the accretion at CE 
phase do not play a very important role on the number of formed 
compact object binaries. 

In models "F" and "L" we consider the results for different 
treatment of the stable MT phases. In particular, in model F2 we 
consider the case of conservative evolution (all mass and angu- 
lar momentum transfered to the companion). Change of the MT 
mode from the non-conservative (standard model) to conservative 
evolution (F2), does not change the numbers by much. The model 
Fl with highly non-conservative evolution, decrease the numbers 
of formed compact object binaries rather significantly, but it is 
rather improbable, since the estimated material loss is probably not 
as high as assumed in model Fl iMe urs & van den Heuvell ll989). 
Since we have adopted quite a large value for specific angular 
momentum of lost material in non-conservative MT episodes, we 
should really concern only the model LI with lowered specific mo- 
mentum leaving the binary. As we see from Table 1 , the numbers in 
model LI are almost the same as for standard model. We thus con- 
clude, that the treatment of non-conservative MT phases does not 
have a great influence on the number of compact object binaries. 

Winds of massive stars may play an important role on the pop- 
ulation of compact objects. In model Gl we decrease all the wind 
mass loss rates by factor of 2. The smaller the winds, the more mas- 
sive compact objects formed, and more BHs are formed as com- 
pared to NSs. However, the total number of compact object binaries 
is basically unchanged in this model. 

The flatter IMF slope of model J increases slightly (as slope 
was not changed by much) the number of heavy stars (progenitors 
of compact objects) and thus leads to a slight increase in the num- 
ber of compact object binaries. The IMF slope for massive stars is 
rather well determined l Kroupa & Weidner 2003), and as expected 
the small change does not affect the population. 

Change of initial mass ratio distribution may have severe ef- 
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feet on the numbers of compact object binaries (models Ml-2). In 
model Ml most of the progenitor systems are formed with ex- 
tremely small mass ratios. Therefore, once they reach first MT 
phase, it is usually dynamically unstable, leading to spiral in and 
merger of components, aborting the formation of compact object 
binary. This explains the small number of formed binaries in model 
Ml, and warrants survival of the systems (large numbers of com- 
pact object binaries) with rather equal masses in model M2. Since 
the initial mass ratio distribution is not easily measurable and con- 
strained, models Ml-2, although rather extreme, should be taken 
account in further analysis. 

In the end, several other models, do not have much influence 
on the production efficiency of compact object binaries. These in- 
clude models with different metallicities (Zl-2), different assump- 
tion on initial eccentricities (S) and finally the model in which we 
change the regime of the fall back in formation of compact object 
binaries (O). 

The masses of the compact objects are strongly affected in 
some of the models. In particular the masses change quite drasti- 
cally in models Gl and G2 where the stellar winds are changed. 
Decreasing the stellar winds (model Gl) allows the star to develop 
more massive cores and consequently leads to higher masses of the 
compact objects formed. A similar effect is connected with decreas- 
ing the metallicity (models Zl and Z2), since lower metallicity stars 
have smaller winds. On the other hand within the model G2, where 
the winds are artificially increased, massive stars do not have the 
time to develop massive cores and no compact objects above 3 Mq 
are produced. The masses of the compact objects are also affected 
by varying parameters in equation Q In model O we increase the 
upper bound of the fall-back range to 14 Mq . this leads to smaller 
masses of the black holes produced as even for high mass cores 
some fraction of the mass is still expelled. The population of com- 
pact object binaries in model El, with reduced common envelope 
efficiency contains additional systems with massive black holes. 
Within model C we turn off the hypercritical accretion onto com- 
pact objects in common envelope events. This primarily influences 
the masses of neutron stars and low mass black holes as they have 
no possibility to increase significantly. 

2.3 Distributions of masses and mass ratios 

In the output we note the masses of the compact objects in each bi- 
nary and the lifetimes: the stellar lifetime from the formation at the 
zero age main sequence to formation of a double compact object, 
and the lifetime as a double compact object binary until it merges 
due to gravitational wave emission. We denote the sum of the two 
lifetimes as the total lifetime of the binary T. 

In this calculation we assume for simplicity that the space is 
Euclidean. We denote the masses of the components in each binary 
m\ and m 2 , the mass ratio is q = m\/m 2 < 1. The formation rate 
of compact object binaries with a given mass ratio q, the mass of 
the primary m 2 (the greater of the two masses), and the lifetime T 
at a given cosmic time t is 



dF{m 2 ,q,T,t) 
drri2dqdT 



S (t ) j 'sin 



(M*) Ntot 

8{m 2 - m 2 )5{q - q)S(T - T l 



(2) 



N CCOB 

E 



where S(t) is the star formation rate at the time t, f S im is the frac- 
tion of stars out of a total population that we simulate, (M t ) is 
the average mass of a binary in the stellar population, and Nccob 



is the number of coalescing compact object binaries formed in a 
simulation of N to t binaries. Our aim is to calculate the observed 
merger rate by an observer on Earth at present, which we denote as 
to- The coalescence rate at a distance r from the observer is then 
given by 



dfcoai(r) _ f , , dF(m 2 , q,t',t -r/c-t') 



dm 2 dq J dm 2 dqdT 
Inserting equation|3|into|2|we obtain 

dfcoaljr) fsim 

dm 2 dq ~ (A/*) 



(3) 



(4) 



NCCOI 

E 



S(m 2 - m 2 )8(q - q')S(t - r/c - T) . 



The observed rate is obtained by integrating equation [4] over the 
volume in which the binaries are observable 



dR 
dm 2 dq 



dV 



dfcoaljr) 

dm 2 dq 



(5) 



IV(m 2 ,q) 

We note that for a constant star formation rate the lifetimes of the 
binaries do not enter the observed rate. 

We first calculate a volume limited distribution of masses, i.e. 
we assume that all binaries coalescing in a given volume V are ob- 
servable regardless of the mass m 2 and the mass ratio q. This cor- 
responds to observing for example the entire population of a given 
galaxy or a galaxy cluster. Here we also assume that the star forma- 
tion history was constant. The normalized volume limited distribu- 
tion of masses and mass ratios is 



P(q,m 2 ) = N C c Ql 



Nccob 



S(q - q')8(m 2 - m 2 ) , 



(6) 



In the case of realistic detectors the volume of integration in 
equation|5|will depend on m 2 and q. Our calculation of the distribu- 
tion of masses and mass ratio relevant for detecting merging bina- 
ries with gravitational waves follows the calculations presented ear- 
lier in l Buli k & Belczvhskil f2003 ) . Here again we assume that the 
star formation rate was flat and that the Universe is Euclidean and 
uniformly filled with stars. The signal to noise in high frequency 
gravitational wave detector s from an inspiral of a stellar mass bi- 
nary is given bv j Fhu^^ ChernofJl993l : lBonazzola & Marckll994L 
iFlanaean & HugheJl99 $T ~ 



(S/N) = ± 



M 



M P 



5/6 



(7) 



where M = {m\m 2 ) 0S {m\ + m 2 )~~ ' 2 is the chirp mass, d is the 
distance, and the Ai depends on the details of a particular detec- 
tor. Thus a coalescence of a binary with a chirp mass M 1 will be 
visible up to a distance proportional to (./VP) 5 / 6 and the volume of 
integration in equation|5|will be V 1 oc (A / ! l ) 5/ ' 2 . The observability 
weighted distribution of masses and mass ratios is therefore 



Pobs{q,m 2 ) = K~ 



' Vi5(q - q')5(m 2 - m 2 ) . 



(8) 



where K = Vj, Again the lifetimes of the binaries do not enter 
the weights in equation [8] because of assumption of constant star 
formation rate history. This distribution is more realistic as it cor- 
responds to a case of an instrument with the sensitivity allowing it 
to detect binary coalescences in a large ensemble of galaxies. Re- 
laxation of the assumption of constant star formation and taking 
into account a realistic cosmological model has been discussed in 
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Volume limited distribution, model A, NSNS binaries 



Flux limited distribution, model A .NSNS binaries 




0.9 1.0 



Volume limited distribution, model C, NSNS binaries 



Flux limited distribution, model C.mmt ,NSNS binaries 




Figure 1. The volume limited and observability weighted distributions of parameters of compact NSNS binaries obtained within model A - top panel and 
model C - bottom panel. The region in dark gray encompasses 68% of the systems, the medium gray corressponds to 95%, and the light gray corresponds to 
all binaries in the simulation. 



iBulik et alliool for the case of the distribution of observed chirp 
masses and was shown not to be significant. 



3 RESULTS 

In the following we will assume that the maximum mass of a neu- 
tron star is 3 Mq . All objects above this value will be considered as 
black holes. In our simulations the minimum mass of a neutron star 
is 1.2 Mq. Thus we can classify all binaries as double neutron star 
(NSNS), black hole neutron star (BHNS), or double black hole bi- 
naries (BHBH). These three categories will be analyzed separately. 
We will present the distributions of binary parameters in the space 
spanned by the mass of the primary (the more massive component 
of a compact object binary) and the mass ratio q. 



3.1 NSNS binaries 

We present the volume limited distributions of q and vn-i in obtained 
in the framework of model A in the top panels of Figure Q The 
volume limited distribution for the case of NSNS systems (top left 
panel) exhibits a peak for the systems with nearly equal masses just 
above the minimal mass of a neutron star. This roughly corresponds 
with the observations of pulsars where most systems have similar 
masses. There is however, a long tail in the distribution extending 
to systems with large mass of the primary, and low mass ratio, these 
are systems consisting of a neutron star with a mass near the max- 
imum value and a companion neutron star with a low mass. Such 
systems have chirp mass about 1.5 times larger than the binaries 
from the above mentioned peak. Therefore in the flux limited dis- 
tribution - top right panel in FigureQ-these low mass ratio systems 
are showing up more prominently. It is however possible that the 
maximum mass of neutron star is lower than 3 Mq and some of the 
systems shown here harbor low mass black holes rather than neu- 
tron stars. However even for the maximum mass of a neutron star 
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Volume limited distribution, model A ,BHNS binaries Flux limited distribution, model A ,BHNS binaries 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Mass ratio Mass ratio 



Volume limited distribution, model O ,BHNS binaries Flux limited distribution, model O ,BHNS binaries 




Tt - 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Mass ratio Mass ratio 

Figure 2. The volume limited and observability weighted distributions of of the parameters of BHNS binaries within model A (top panel), model O (middle 
panel) and model Z2 (bottom panel). The region in dark gray encompasses 68% of the systems, the medium gray corressponds to 95%, and the light gray 
corresponds to all binaries in the simulation. 
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Volume limited distribution, model A ,BHBH binaries 



Flux limited distribution, model A ,BHBH binaries 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Mass ratio 




Volume limited distribution, model O ,BHBH binaries 



Flux limited distribution, model O ,BHBH binaries 





0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Volume limited distribution, model E1 ,BHBH binaries 



Flux limited distribution, model E1 ,BHBH binaries 





0.3 0.4 0.5 0.6 0.7 0.8 

Mass ratio 



0.9 1.0 




Figure 3. The volume limited and observability weighted distributions of of the parameters of BHBH binaries within model A (top panel), and model O 
(bottom panel). The region in dark gray encompasses 68% of the systems, the medium gray corressponds to 95%, and the light gray corresponds to all binaries 
in the simulation. 
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of 2.0 Mq there is still a large fraction of low q systems in the flux 
limited distribution. We have examined all 21 model listed in Table 
1 and nearly all the models show a similar pattern in both the vol- 
ume limited and flux limited distributions. There is one exception - 
model C - for which we present the relevant distributions in the bot- 
tom panels of FigureQ In this model we turn off the possibility of 
hypercritical accretion onto a compact object in common envelope 
phase. This effectively shuts off the possibility of increasing signif- 
icantly the mass of a neutron star through accretion. Therefore in 
the volume limited distribution there is quite a large concentration 
of systems with both masses below 1.5 — 1.6 Mq , and little number 
of binaries with low mass ratios. Consequently in the flux limited 
distribution the binaries with low mass ratio are nearly absent in 
contrast to model A. However, the 68% contour includes systems 
with q > 0.75 and 7772 ~ 2.0 Mq, as well as some binaries with 
q ~ 0.75 and the mass of the primary close to the maximal mass 
of a neutron star in our model. 



3.2 BHNS binaries 

We present the volume limited and the flux limited distributions 
of q and 777,2 in Figure [2] The top panel of Figure [2] corresponds 
to the standard model A. The volume limited distribution shows a 
large number of binaries along a stripe stretching from q w 0.7 and 
m2 ~ 4 Mq to q ~ 0.2 and 7712 ~ 10 Mq. Systems above and to 
the right of this stripe would be classified as BHBH binaries. The 
volume limited distribution is dominated by binaries with low mass 
black holes. In the flux limited distribution the binaries with higher 
mass black hoes and low mass ratios start to play an important role. 
This is due to the balance between the falling mass function and the 
increase in the chirp mass with increase of the mass of the primary. 
For nearly all models this leads to dominance of binaries where the 
black hole primary has a mass between 6 and 12 Mq, and the mass 
ratio is somewhere from 0.3 for the most massive black holes to 
0.5 for the moderate mass ones. 

After examining the 21 models of Table 1 in the case of BHNS 
binaries one can distinguish two other classes of models with dif- 
ferent distributions of binary parameters. The first class consists of 
models , L2, M2, and O. For this models we show a representa- 
tive case (model O) in the middle panel of Figure [2] The common 
characteristic of this class of models is that the population of black 
holes in binaries lacks the very massive ones for various reasons. 
In the case of the model O shown in Figure |2| this is because we 
increase the range of masses of fall-back formation of black holes. 
Within the models El, L2, M2 such binaries have a smaller chance 
of formation because of altering the treatment of the mass transfer 
events, in consequence the is very little low mass ration BHNS bi- 
naries and the volume limited distribution is very similar to the flux 
limited one. 

A separate class of consists of models El, Z2 and Gl. In their 
case high mass black holes are easily formed, because of decrease 
of strength of stellar winds (Z2 and Gl). Model El favors sur- 
vival of systems with high mass first born stars and leads to ef- 
fective production of extreme mass ratio compact object binaries. 
We present the two distributions for the case of model Z2 in the 
bottom panel of Figure |2| Here the volume limited distribution is 
dominated by a nearly vertical stripe at q ~ 0.2. In the flux limited 
distribution the binaries with high mass primaries 777,2 ~ 12 Mq 
and 0.1 < q < 0.3 are dominant simply because of their high 
chirp mass. 



3.3 BHBH binaries 

The case of BHBH binaries is presented in Figure|3| The top panel 
corresponds to the model A. The volume limited distribution fills 
more or less uniformly the region allowed for the BHBH binaries. 
In the flux limited distribution there is a preference for the high 
mass ratio (nearly equal mass) and high mass ratio systems, i.e. 
these filling the top right corner of the plot. Thus the flux limited 
distribution is dominated by the systems with q < 0.6 and 7772 near 
the maximum mass produced in a given model. All models seem to 
follow this general trend. 

In models O the maximum mass of a black hole is decreased. 
We present the results of a calculation using model O, in the middle 
panel of Figure[3] In these case the volume limited distribution also 
fills nearly uniformly the region allowed for black holes. However 
because of it smaller size the range of chirp masses for given mass 
ratio is not as large and the flux limited distribution is only slightly 
shifted to higher masses in respect to the volume limited one. 

Another special case - model El - is presented in the bottom 
panel of Figure [5] Here because of lowered CE efficiency forma- 
tion of equal mass compact object binaries is favored, while ex- 
treme mass binaries are preferentially formed. The volume limited 
distribution in this case is dominated by systems with mass ratio in 
the range 0.4 < q < 0.6 and a tail extending to q = 0.8. 

Models of the class that favors production of massive black 
holes (Z2, Gl) do not lead to qualitatively different results than 
model A. In these models the flux limited distributions tend to con- 
centrate around binaries with higher mass ratio, q > 0.7 and higher 
total masses than in the standard model A. 



4 CONCLUSIONS 

We have calculated the expected distributions of masses and mass 
ratios of compact object binaries to be observed in gravitational 
waves. The results are based on the Star Track binary population 
synthesis code. For most of the models the observability weighted 
distribution of double neutron star systems has two peaks: one with 
the mass ratio almost unity and both masses near the smallest mass 
allowed for neutron stars, and another with small mass ratio, con- 
sisting of stars with the mass near the maximum mass of a neutron 
star in a binary with a star close to the minimum mass. The real- 
ity of this second peak depends on the assumed maximum mass of 
a neutron star: the lower the maximum mass of a neutron star the 
smaller the small mass ratio peak. The distribution of black hole 
neutron star binaries peaks at mass ratios between 0.3 and 0.5. The 
bulk of observed double black hole binaries has mass ratios above 
0.7. We have shown that these results are rather generic and depend 
weakly on the choice of a particular model of stellar evolution. 

The crucial parameter determining the shape of the distribu- 
tion of the observed NSNS binaries is inclusion of the hypercritical 
accretion onto compact objects in common envelope events. In the 
case of BHNS and BHBH binaries the most important parameters 
are these that alter the masses of the black holes in such binaries, 
and the common envelope efficiency. The masses may be altered 
either due to the mechanism of compact object formation in su- 
pernova explosions, or due to particular treatment of mass transfer 
events. The distribution of the BHNS binary parameters is most 
sensitive to these changes. However, we must note that hat the ob- 
served sample is dominated by the BHBH binaries. For most mod- 
els mo re than 90% of observed s ystems are double black hole bi- 
naries <Bulik & Belczvnskil2003l) . 
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These results can be used as a guideline for choosing the ini- 
tial conditions in numerical simulations of mergers of compact ob- 
ject binaries. Additionally the results of this work can be used in 
preparing data analysis software using templates for detection of 
gravitational waves from compact object inspiral. Coalescences of 
BHBH binaries dominate the observed sample, and we find that the 
observability weighted distribution is peaked around nearly equal 
mass binaries. We find that in most models the flux limited sam- 
ple of NSNS binaries contains a large fraction of non equal mass 
objects. We conclude that the initial search for gravitational waves 
from coalescences of compact object binaries should concentrate 
on BHBH coalescences with mass ratio close to unity, and the low 
mass ratio NSNS coalescences should be taken into account. 

Finally, we note that this work only includes binaries that 
evolved in galaxies, and neglects all possible effects, like multiple 
stellar interactions that are relevant for evolution in dense stellar 
clusters. 
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